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Abstract. Quantum computation, in particular Grover’s algorithm, has aroused a 
great deal of interest since it allows for a quadratic speedup to be obtained in search 
procedures. Classical search procedures for an N element database require at most 
O(N) time complexity. Grover’s algorithm is able to find a solution with high prob¬ 
ability in 0(y/ r N') time through an amplitude amplification scheme. In this work we 
draw elements from both classical and quantum computation to develop an alternative 
search proposal based on quantum entanglement detection schemes. In 2002, Horodecki 
and Ekert proposed an efficient method for direct detection of quantum entanglement. 
Our proposition to quantum search combines quantum entanglement detection along¬ 
side entanglement inducing operators. Grover’s quantum search relies on measuring a 
quantum superposition after having applied a unitary evolution. We deviate from the 
standard method by focusing on fine-tuning a unitary operator in order to infer the 
solution with certainty. Our proposal sacrifices space for speed and depends on the 
mathematical properties of linear positive maps A which have not been operationally 
characterized. Whether such a A can be easily determined remains an open question. 
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1 Introduction 

Computer scientists are often faced with the task of constructing algorithms capable of deliv¬ 
ering a solution for a given problem. For some problems it is possible to engineer algorithms 
capable of producing a solution with a number of computational steps that is bounded by a 
polynomial n k where n is the length of the input and k some constant. The class of problems 
for which a polynomial-time algorithm exists is known as P. Problems belonging to P are 
usually seen as being efficiently solvable, i.e. tractable. Class EQP represents the quantum 
equivalent of P. 

For other problems it is possible to verify in polynomial-time if a given configuration is a 
solution, although there are no known methods for efficiently calculating a solution. For these 
type of problems, there is no alternative but to perform an exhaustive search of all possible 
configurations. The class NP consists of those problems whose possible configurations can be 
verified in polynomial-time. Clearly, P C NP since the possibility of constructing a solution in 
polynomial time also implies that a solution can be verified efficiently. One of the outstanding 
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questions in computer science consists in determining if the class NP is equivalent to the class 
P, i.e. P=NP? Traditionally, approaches to answering this question have focused at a subclass 
of NP, namely NP-complcte problems. This subclass contains those problems which are both 
NP and NP-hard. A problem is said to be NP-hard if an algorithm capable of solving it can be 
translated into an adequate algorithm for any NP problem. By its own definition, an efficient 
solution for a problem in NP-complete implies that an efficient solution exists for all problems 
in NP. 

The first clues that some problems which are classically hard may have an efficient quantum 
solution were provided in T|. Shor’s algorithm for efficient factorization [2] reinforced this 
idea. Later, Grover’s search algorithm [3] provided an asymptotical quadratic speedup over 
classical strategies. The quantum search algorithm systematically increases the probability of 
obtaining a solution with each iteration. After the algorithm has concluded, a measurement 
is performed in a quantum superposition, in order to obtain a solution with high probability. 
The superposition state represents the set of all possible results. Grover’s approach sparked 
interest by the scientific community on whether it would be possible to devise a faster search 
algorithm. Unfortunately, it was proved that the search problem cannot be solved under 
{2('/N") time [4] using standard quantum computation approaches. 

In this work we present an alternative search method based on the principles of tree search de¬ 
composition and quantum entanglement detection. Unlike traditional approaches, we opt not 
to concentrate our efforts on measuring a quantum superposition of possible values. Rather, 
we are more interested in exploiting the unitary operator that is applied to a quantum super¬ 
position in order to infer possible solutions with certainty. However, an implicit caveat exists 
associated to our quantum search proposal. Namely, our system implies a tradeoff between 
speed and space that will become apparent in the following sections. 

The next sections are organized as follows: Section[2]focuses on presenting the details of an NP- 
complete problem, namely the Boolean satisfiability problem, alongside classical tree search 
techniques of examining the problem space. Section[3]presents our hybrid approach, combining 
tree search decomposition alongside with quantum entanglement detection schemes. Section 
|4] presents the conclusions of this work. 


2 Traditional approaches to tackling NP-Complete problems 

The satisfiability (SAT) problem was the first problem ever shown to be NP-complete [S]. SAT 
asks whether a given boolean formula is satisfiable. Any polynomial-time algorithm capable 
of solving SAT automatically enables an efficient solution for all of NP. In complexity theory, 
the satisfiability problem is a boolean formula <j> composed of [Bj 

— n boolean variables: aq, aq, • • • ,x n ] 

— m boolean connectives: any boolean function with one or two inputs and one output, such 
as A (AND), V (OR), -i (NOT), —> (implication) and e->-(if and only if); 

— parentheses. 

We are interested in determining a set of values for the variables of (/>, i.e. variable configura¬ 
tion, which cause the overall expression to be satisfiable, i.e. evaluate to true. At any given 
point in time we need to consider the n variables alongside m gates, i.e. we can verify any 
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configuration in n + m time. However, the number of possible configurations to consider grows 
exponentially with the cardinality of the variable set. As an example lets consider the simple 
formula presented in Expression [T] 


4> = (%i A x 2 ) V x 3 (1) 

The standard approach to solve such a problem would be to enumerate all possible configu¬ 
rations of the m variables. This procedure can be better understood with the help of a simple 
tree diagram such as the one illustrated in Figure [Q At each depth level a specific the pos¬ 
sible values for a specific binary variable are considered, e.g. depth 0 considers the possible 
values for X\. depth 1 considers variable x 2 and so on. With each depth level an additional 
binary variable is taken into account. Considering n binary variables requires examination of 
2" possible leaf states, i.e. !?( 2 n ). In tree search vocabulary these states are also known as 
paths. If the specific case of Expression |T] is mapped into the tree elements of Figure [T| then 
paths 1, 4, 6 , 7 and 8 would evaluate to true. 



3 Approach 

How can we proceed by developing an alternative approach to that of Grover’s? First lets 
start by considering the following scenario: suppose we have a bipartite quantum system 
respectively labeled as the query register, |g), and the answer register, |a), acting on Hilbert 
space TL = H q ® 'Hq. The query register is an n-qubit register where possible values for the 
binary variables of the SAT problem will be setup, i.e. | q) = \x\x 2 • • ■ x n }. Notice that in order 
to gain a quantum advantage over classical computation we need to place |g) in a uniform 
superposition of the computation basis. This can be done efficiently by applying the Hadamard 
transform H a total of n times to the n-qubit state |0), i.e. H® n |0}®" = - 7 =^ \ x )• Such a 

procedure enables the creation of a superposition containing an exponential number of states, 
each of which representing a possible tree path, by only employing a polynomial number of 
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gates. The answer register contains a single qubit which is initialized to state |0). The overall 
state of the system can thus be described as illustrated in Expression [2] 


k)|o) 



2 n — 1 


E i*>i°> 

x=0 


( 2 ) 


Additionally, suppose that a quantum oracle with the form presented in Expression [3] is con¬ 
structed. The auxiliary function </>: {0,1}™ —>-{0,1} employed simply verifies if an argument 
is a solution or not for a specific SAT instance, as illustrated by Expression 2J We should be 
careful to point out that an efficient oracle responsible for verifying the validity of a variable 
configuration for a specific <f> can be easily constructed by mapping the m boolean connectives 
of the network onto a reversible circuit (see (3) in order to ensure a unitary mapping. 


O |g)|a) = |q)|a ® </>(<?)} 


( 3 ) 


4>(q) = <j>(x>x n ) 


1 if q evaluates to true 
0 otherwise 


( 4 ) 


If oracle O is applied to the combined state of Expression [2] a result like the one illustrated in 
Expression [5] may be obtained, where \i/j') denotes the overall superposition evaluation. For 
simplification issues we assume that there exists at least a solution. Naturally, some of the 
query values produce a solution, whilst others do not. 


W) 


' | 00 • • • 0 )| 0 > 


bits 



2 n — 1 


E O|*)|0) 

21—0 


= 


|00 • • • 1 )| 0 ) 
111 ... 0 ) 11 } 


1 | 11 -- 1 )| 0 ) 


( 5 ) 


From this point on the system’s state can no longer be expressed as a tensor product be¬ 
tween query and answer register, i.e. the system becomes entangled. Quantum entanglement 
is a key feature of quantum mechanics which details the connections between subsystems of 
compound quantum systems. It was a key aspect of the quantum world formalism proposed 
by von Neumann in 1932 0. Although the intriguing impacts of quantum inseparability were 
only later grasped by Einstein, Podolsky and Rosen 0 alongside Schrodinger m- Quantum 
entanglement is also a key resource in quantum information. 

Mathematically, we can describe the state of each register by tracing out the remaining reg¬ 
ister, through the partial trace mechanism. In this case we are interested in the overall state 
of the answer register. In order to calculate the partial trace of the answer register we first 
need to calculate g the density operator of the quantum state presented in Expression [5] The 
overall form for g a is illustrated in Expression 0 
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Q= \lf)){lj}\ 

= 7=r(|00---0>|0} + |00---l>|0> + - 
-^((00 • • • 0 |( 0 | + (00 • • • 1 |( 0 | + • 
= Jr|00---0)|0)«00---0|(0| + (00-- 
Jr | 00 -- ■ 1 )| 0 >(< 00 - ■ - 01(01 + ( 00 -• 
+ • • • ~f" 

^ r|ll - - • 0 )| 1>(<00 • • • 0 |< 0 | + <00 - • 
^ r|ll - - • 1 }| 0>(<00 • • - 0 |< 0 | + <00 - • 


- + 111 -- 

• 0 )| 0 } + 

111 ... 

1 >| 0 » 


• + < 11 " 

■- 01(01 + 

( 11 ... 

1 |< 0 |) 


1 |< 0 | + - 

•■ + < 11 " 

- 0 |( 0 | 

+ ( 11 " 

- 1 |< 0 |)+ 

1 |( 0 | + - 

•■ + < 11 " 

• - 0 |( 0 | 

+ ( 11 " 

- 1 |< 0 |)+ 

1 |( 0 | + - 

•■ + ( 11 " 

• - 0 |( 0 | 

+ ( 11 " 

- 1 |< 0 |)+ 

1 |< 0 | + - 

•■ + ( 11 " 

• - 0 |( 0 | 

+ ( 11 " 

- 1 |< 0 |)+ 


( 6 ) 


Q a = Tr,(e) 

= £((00 - - - 0||00 - - - 0 }| 0 )( 0 | + <00 • • • 1||00 • • • 1 )| 0 )( 0 | + 

H-b (T) 

(11 • • • 0||11 • • • 0 )| 1 }( 1 | + (11 • • • 1||11 • • • 1 )| 0 )( 0 |) 

= £[(2"-l)|0)(0| + |l)(l|] 


Generally, the result presented in Expression [G] can be improved if we take into account the 
number of solutions. Accordingly, let k denote the overall number of solutions, then g a takes 
the form shown in Expression^ Notice that the overall state is separable only when k = 0, i.e. 
no solution exists, or when k = 2 n , each value belonging to [0, 2 n — 1] is a solution. Otherwise, 
the system is entangled. 


^ = i r [(2"-fc)|0)(0|+fc|l)(l|] (8) 

Thus, the problem of determining whether or not a solution to a problem exists can be 
reduced to the problem of determining whether the overall quantum state is separable or 
entangled. 


3.1 Quantum entanglement detection 

The quantum separability problem consists in determining if a given a density matrix g rep¬ 
resenting a quantum state is entangled or separable m Efficiently deciding on the nature of 
such states has grabbed researchers attention and remains a problem of crucial importance to 
the fields of quantum computation and information m ■ Generally speaking, quantum entan¬ 
glement is studied in accordance with a varied mix of properties (just to name a few of these: 
bipartite vs. multipartite systems, pure vs. mixed states, bound entanglement; for exhaustive 
reviews please refer to na, m and m- It is important to mention that the quantum sep¬ 
arability question has been approached from the classical and quantum perspectives. These 
approaches typically consider the nature of the input (classical vs. quantum), and whether 
any required processing will be performed on a classical or quantum computer m • This prob¬ 
lem was shown to be NP-hard classically m- However, as mentioned in m the processes 
involving both quantum input and processing have not been thoroughly investigated. 

In the case of our specific approach we would only need to consider bipartite quantum systems 
with mixed states. As pointed out in |l8i the mixed state requirement stems from the fact 
that any potential laboratory demonstration of this approach would have to deal with mixed 
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states rather than pure ones, due to the uncontrolled interactions with the environment. 
These requirements are present in one of the existing quantum detection schemes, namely 
the one proposed in [19j . The method employed by the authors is experimentally viable and 
provides for a direct detection mechanism of quantum entanglement. Their approach is based 
on the theoretical foundations laid down in [18] . The method determines whether a state g 
is separable or not, i.e. entangled, based on the mathematical properties of linear positive 
maps acting on matrices. More specifically na, let Md —> Md be the space of matrices of 
dimension d, a map A : Md —> Md is called positive if it is Hermitian and has non-negative 
spectrum. Additionally, the map A is completely positive if and only if I ig) A is positive for 
identity map I on any finite-dimensional system. A state g is separable if and only if the 
result presented in Expression [5] is observed for all positive but not completely positive maps 

A : Md —> Md- 


[/®A](p)>0 (9) 

Expression [S] cannot be directly used since it requires knowing state g beforehand. Addition¬ 
ally, positive maps A cannot be directly implemented in laboratory. Fortunately, it is possible 
to obtain a physically realizable map by mixing an appropriate proportion of [I ® A] with a 
depolarizing map. This approach allows for a new map [I ® A] to be obtained, which have 
been referred to as structural physical approximations. For more on this subject please refer 
to [5^. The separability criterion can then be restated as follows m : g is separable if and 
only if for all positive maps A the condition presented in Expression 1101 is observed. 

^ do) 

Where A corresponds to the most negative eigenvalue obtained when the induced map [(/ (§) 
/) ® (/ ® A)] acts on the maximally entangled state of the form A? 1 l*)l*)- Accordingly, 
Expression [10] states that the lowest eigenvalue of the transformed state g' — [I ® A](g) should 

j2 \ 

be greater than d 4 X ^_ 1 for g to be separable. 

The authors devised a method which allows for an estimate of the lowest eigenvalue to be 
obtained efficiently and directly. It requires that a joint measurement be performed on N 
copies of state g'. The overall input density operator of the estimation problem is g'® N , which 
exists on the TVth tensor power [2T|. The error e associated with the estimate of the 

lowest eigenvalue decreases exponentially with N. Such a measurement can be represented 
as quantum network implementing projections on the symmetric and partially symmetric 
subspaces ESI • An efficient method addressing these questions was proposed in m requiring 
a number of auxiliary gates that grows quadratically with the dimension of the input, i.e. 
0(n 2 ), where n is the number of bits. If g' represents the state of an n qubit register, then each 
additional tensor power will mean that another n bits should be accounted. Consequently, 
an g'® N system will have a total of N x n bits. Which means that the quantum network 
responsible for estimating the lowest eigenvalue will have 0{N 2 n 2 ) complexity. 

Clearly, this approach is dependent on map A which have not been operationally characterized 
so far [23]. As pointed out in (14] in general the set of positive but not completely positive maps 
is not characterized and it involves a hard problem in contemporary linear algebra. However, 




Can Quantum Entanglement Detection Schemes Improve Search? 


7 


for low dimensional systems, namely those with dimension 2 ® 2 or 2 ® 3, the positive partial 
transpose map proposed in m can be employed as the A. In m the authors draw attention 
to the fact that ‘Recently, the progress in this direction has been made [25] [26] which suggests 
that tests of separability based on positive maps will soon acquire practical meaning beyond 
the scope of two-qubit systems.” Whether such a map A acting on 'Hd^'Hd quantum systems 
can be determined remains an open question. 


3.2 Subset entanglement inducing oracle 


The fl(VF) lower bound for quantum search employing oracles working on the full range of 
searchable items implies that an alternative search approach has to be devised. In classical 
tree search it is a standard technique to start by analyzing subtrees and deciding whether 
these may eventually lead to a solution. Based on problem requirements it is possible to 
automatically exclude, i.e. prune, certain subtrees. The act of pruning may eventually be 
responsible for large sections of the tree to be discarded, and therefore allow the search to 
terminate faster. We will draw inspiration from these concepts of classical search in order to 
develop our approach to quantum hierarchical search. 

Quantum algorithms employing traditional oracles provide at most a polynomial advantage 
over classical algorithms for total functions, i.e. functions defined for the whole of { 0 , l} n , 
where n is the number of bits. The oracle model contemplates syperpolynomial advantage 
but only when partial functions are defined which operate on a subset of {0,1}" 22 • No¬ 
tice that classical search can be viewed as a procedure which evaluates subsets of an initial 
range. Since in quantum computation the oracle operator can be applied to a superposition of 
computational basis, evaluating subsets is equivalent to only evaluating specific ranges of the 
superposition. Accordingly, it is possible to develop an oracle responsible for evaluating only 
a certain subset of the initial range [0,2" — 1] allowed with n qubits. Although we are only 
interested in evaluating specific subset there are other alternatives for trying to decompose 
a quantum search space. For instance, Grover concluded in [281 that determining the first n 
bits of a solution by employing amplitude amplification schemes is only slightly easier than 
determining the total bits. 

This model for a range specific entanglement inducing oracle can be described as presented in 
Expression IIII which employs an auxiliary function /[ a ,&](<?) defined in Expression 1121 In the 
case of the SAT problem it would be convenient to define /[ a ,&](<?) as 4>\ a ,b\ (<?)• 

0[ a ,b\\q)\a) = k)|a© /[a, {,](<?)) (11) 


f[a,b] (?) 


1 if f(q) is a solution and q £ [a, b\ 
0 otherwise 


( 12 ) 


As was previously pointed out, the oracle evaluation process has the overall effect of making 
the quantum registers entangled. By testing whether the oracle has induced, or not, quantum 
entanglement it is possible to decide if the search procedure should continue decomposing a 
range, or if another range should be tested. Ideally, the entanglement detection scheme should 


L. Tarrataca, A. Wichert 


present some type of polynomial upper-bound behavior such as the one described in the pre¬ 
vious section. □ If the state g resulting from applying an oracle O with the form presented 
in Expression [12] is separable then the range evaluated can automatically be discarded. Dis¬ 
carding a wide range of potential candidates en masse can be understood as the classical tree 
search operation of pruning certain subtrees. On the other hand, if g is entangled then it is 
possible to further decompose the associated range. Eventually, this sort of recursive branch 
and bound procedure, by constantly readjusting the range of oracle O, will “zoom in” on 
a solution. Extending this method to problems presenting multiple solutions would require 
systematically focusing on previously non-expanded but solution-containing ranges. 

Notice that this approach requires a new oracle to be defined with each iteration in terms of 
a specific subset that may be entangled. The set of oracles applied throughout the search can 
be viewed as a single “dynamic” oracle, which differs substantially from the standard “static” 
oracle applied in quantum search. Additionally, in contrast with Grover’s algorithm, we are 
not interested in performing an amplitude amplification process, but rather we are concerned 
with decomposing the quantum search space. 


3.3 On the growth of the number of copies required 

Clearly, one question still lingers: What can be said about the number of copies N of the 
system that are required? According to m any procedure that on input | ipz) guesses whether 
Z = X or Z = Y will guess correctly with probbability at most 1 — e = | + \\Jl — 6 2 ' where 
6 = | (ipx 11 V’y) I ■ For our particular case we are interested in distinguishing two specific cases, 
namely: 

— \4>solution) which results from applying U\ijj) when one solution exists; 

— no—solution) which results from applying U\ip) when no solution exists; 

For search spaces of dimension L the initial amplitudes at associated with each computational 
basis i of the superposition \ip) is After having applied oracle O the two states remain 
exactly equal except for two computational basis where the amplitudes permuted. This means 
that when calculating the inner product the permuted computational basis will sum up to 
zero. Accordingly, the inner product will sum the value —^ a total of L — 1 times, i.e. 

& — (^solution | |^no— solution) — ^ (13) 

Given a tensor product of N items, Expression [T3] evolves into Expression [TT] The three- 
dimensional plof of 5® N as a function of L € [2 1 ,2 30 ] and N £ [2 1 ,2 30 ] is illustrated in Figure 

El 


T _ 1 N 

S® N = {^solutionUno-solution)® N = (14) 


* The entanglement detection approach described in [T9] requires the overall bipartite system to be 
d ® d. Consequently, the answer register |a) should have the same dimension than |g), i.e. n bits. 
This requirement has no direct consequences in the overall oracle since unitary evolution can still 
be assured. 
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In order for these states to be distinguished with significant probability the inner product 
5® N presented in Expression 1141 must be made small. However, in order to achieve this one 
needs to chose a number of copies N that grows in accordance with the dimension of the 
search space L , i.e. N = O(L). Consequently, this approach would not provide for any gains 
over classical search. 


3.4 Consequences for efficient entanglement detection schemes 

What would be the consequences if the number of system copies N was not a function of 
the search space? Suppose the proposed search procedure is executed on n-qubits placed on 
a superposition. Initially, the algorithm has to decompose the [0, 2" — 1] initial range. Lets 
assume that any specific range being considered is split in half. Accordingly, the procedure 
needs to verify if evaluating the elements in [0,2 n_1 — 1] produces an entangled quantum 
state g. If this is found to be true then subset [0,2" _1 — 1] can be also split in half and 
evaluated. Otherwise, subset [2" _1 — 1,2" — 1] needs to be decomposed. Independently of 
what subset induces entanglement, the algorithm is able to prune half of the 2" initial states, 
i.e. 2"/2. Accordingly, for iteration i, the oracle is able to focus on 2"/2* states. Clearly, 
when i = n a single state is being considered and consequently a solution can be determined 
with certainty by employing 0(n ) oracle queries. Associated with each oracle query is the 
quantum entanglement detection scheme bringing the overall complexity of our approach to 
0(N 2 n 3 ). In the case of the SAT problem we have to consider the costs associated with each 
oracle query, respectively n + m. Consequently, a solution for SAT would be calculated in 
0(N 2 n 4 + N 2 n 3 m) quantum polynomial time. 

It is our believe that it is not possible to efficiently detect quantum entanglement. If we take 
into account the simplicity of the search procedure designed in Section 13.21 then if such a 
method existed we could efficiently search, i.e. in quantum polynomial time, a problem space 
of dimension d. Accordingly, we can define the following conjecture. 

Quantum entanglement detection conjecture - It is not possible to efficiently detect 
quantum entanglement non-classically since this would automatically imply that a simple 
algorithm exists proving that NP=EQP. 
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The above conjecture stresses the notion that there appears to be a relationship between 
entanglement detection and search in terms of computational complexity. 


4 Conclusions 

The general characterization of positive but not completely positive linear maps A alongside 
quantum entanglement detection schemes and partial range entanglement inducing operators 
may eventually be responsible for producing efficient algorithmic solutions capable of search¬ 
ing exponential-growth search spaces. Although some research has already been carried out, 
further thorough analysis into the subject is still required. However, given that N = 0{L) 
current methods cannot be employed in order to speed up quantum search. 
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